Detection of large acoustic energy flux in the solar atmosphere 
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Tj- ■ We study the energy flux carried by acoustic waves excited by convective 

Q\ • motions at sub-photospheric levels. The analysis of high-resolution spectropo- 

larimetric data taken with IMaX/SlJNRlSE provides a total energy flux of ~ 6400- 
7700 Wm -2 at a height of ~ 250 km in the 5.2-10 mHz range, i.e. at least twice 
the largest energy flux found in previous works. Our estimate lies within a factor 
of 2 of the e nergy flux needed to balan ce radiative losses from the chromosphere 



Oh- 



O 



according to lAnderson fc Athayl (jl989l ) and revives interest in acoustic waves for 



transporting energy to the chromosphere. The acoustic flux is mainly found in 
the intergranular lanes but also in small rapidly-evolving granules and at the 
bright borders, forming dark dots and lanes of splitting granules. 

Subject headings: Sun: photosphere — Sun: chromosphere — Sun: oscillations 
- techniques: high angular resolution — techniques: spectroscopic 



1 The National Center for Atmospheric Research is sponsored by the National Science Foundation. 



-2- 



Introduction 



The long-running debate on the heating of chromospheric layers by acoustic waves is 
fed by the fact that even the largest measured energy fluxes (found by Bello Gonzalez et 
al. 2009,2010) are nearly an order of magnitude s maller than the amount n eeded to balance 
the chromospheric energy losses of 14000Wm -2 (lAnderson &: Athayill989l ). Bello Gonzalez 
et al. found an acoustic energy flux of ~3 000 W m -2 at 250 km and of ~2 000 W m~ 2 
at 500 km from velocity fluctuations measured from narrow-band (FWHM = 1.8 pm) data 
at ~0'/4 spatial resolution on the Fel lines at 557.6 nm and 543.4 nm, formed in mid- 
photospheric and low-chromosphe r ic laye rs, re spectively. Thi s is a f actor of 4-6 larger than 
that given by iFossum fc Carlssonl (120061 ) and ICarlsson et al.l ( 120071 ) from intensity fluctua- 
tions in continuum bands at 160 nm (of ~ 1" spatial resolution) and in CallH formed at 
height s of 430 km and 200 km, respectively. We refer the reader to iBello Gonzalez et al. 



( 120091 ) for a detailed introduction to the acoustic-wave heating debate. 

In this contribution, we present results on acoustic waves from data taken in the Fel 
line at 525.02 nm with the highest spatial resolution considered so far. 



2. Observations and data analysis 



The data an alysed in this study were taken on June 9, 2009 with the IMaX 2D- 



spetr opolarimeter ( IMartmez Pillet et al.ll2010l ) onboard Sunrise (IBarthol et al.ll2010l ; ISolanki et al 
20101 ). This study is performed on a time series of ~23min with 33 s cadence, taken 
from 00:36 UT on in full-polarimetric mode, at quiet disc centre. The data consist of nar- 
row (FWHM = 6 pm) spectropolarimetric filtergrams taken at [±8, ±4, +22.7] pm around 
the Fel 525.02 nm line (Lande factor g = 3). The capabilities of ret rieving the physical 



param eters with the given spectral sampling have been studied by lOrozco Suarez et al. 



( 120101 ) . The data were calibrated according to IMartmez Pillet et al.l (120101 ) (Sect. 9.2) and 
possess a spatial r esolution of 0^15-0^18 , achieved by combini ng; the onboard image sta- 
bilizati on system dGandorfer et al.l l2010i : iBerkefeld et al.l |2010[ ) and phase-diversity tech- 
niques (IMartmez Pillet et al.ll2010l ). The field of view (FOV) considered is 44.8x44.8 arcsec 2 
(815x815 pix 2 ). The line-of-sight (LOS) velocities were determined by measuring Doppler 
shifts of the min imum position of a Gaussian fitting the 4 spectral points within the Fe I line 
plus continuum (IMartmez Pillet et al.ll2010l ). 
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2.1. Response functions 



We base our study of wave phenomena on the analysis of fluctuations in the measured 
LOS velocities. For a better comprehension of the atmospheric layer where the Fel line 
at 525.02 nm is sensitive to velocity fluct uations, we calcu lated veloci ty response functions 
(RF V ) , applying the method described by Elbe et al.l (120011 ) . Following iBello Gonzalez et al. 
( 120101 ). we considered granule (GR) and in tergranule (IGR) model atmospheres based on 
3D simulations from lAsplund et al.l (120001 ) . Results are shown in Fig.[]Ji. In the case 
of a spectrometer with infinite spectral resolution, the RF V calculated for the line min- 
im um in both atmospheres are centre d around 400 km height, in agreement with results 
by IShchukina fc Trujillo Buenol (120011 ). However, the RF V are shifted to lower layers and 
partly broadened when (a) the spectral transmission of the IMaX instrument and (b) the 
Gaussian-fit used in the determination of the line-minimum velocities, are applied to the 
emergent intensities in order to mimic the measurements from IMaX observations. For GR, 
the RF V is centred at ~ 325 km and for IGR at ~ 175 km. 



2.2. Atmospheric transmission to velocity amplitudes 

The broad extent in height of the RF V indicates the large atmospheric range where 
the measured velocities are referred to, i.e. indicates the difficulties in detecting the signal 
of small-scale velocity fluctuations along the LOS. It is then advantageous to estimate the 
transmission of the solar atmosphere to (wave) velocity amplitudes, i.e. the ratio of the 
observed to actual velocities, taking into account the spectral characteristics of our optical 
system and the velocity determination method. 

FigU 

The transmission function is modeled by computing velocities from synthesised IMaX 
profiles from a model atmosphere including velocity fluctuations (waves) along the LOS. We 
consider both the GR and IGR atmospheres following two approaches (i) waves with height- 
independent amplitude and (ii) with exponentially increasing velocity amplitudes. The lat- 
ter i _jmmicj_an_a^ou^Hc eneiw flux constant through the stratified atmosphere. We refer 
to IBello Gonzalez et al.l (120091 120101 ) for detailed descriptions on the transmission function 
calculations and the determination of group velocities derived from the dispersion relation 
of waves. 



Figure[Tb depicts the resulting functions. They show similar transmissions for the GR 
atmosphere when considering height-independent or exponentially increasing velocity ampli- 
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tude. The IGR atmosphere shows a lower transmission for waves with height-independent 
amplitude. Lower transmission means larger correction in the determination of the energy 
flux (see Sect. 13. ip . 



3. Results 

We performed Fourier and wavelet analyses of the velocity fluctuations over the time 
series, as in Bello Gonzalez et al. (2009, 2010). The temporal power spectra were calculated 
at each pixel and averaged over the FOV. A separate study on granular and intergranular 
areas has been performed. Granules (GR) are considered to be those regions in the continuum 
intensity maps where /coat /(/rant) > 1 an d intergranules (IGR) where leant/ (leant) < 1- 
Granular dark dots and forming lanes are thus treated as intergranules. 



3.1. Acoustic energy flux 

The total flux of acoustic energy (per element of area on the Sun) is estimated from the 
spatially averaged Fourier power P v (vi) by 

i r ac,tot(fi) = P P v{Vi) ■ V &ouv {Vi)/T(Vi) ■ AUi , (1) 

i 

with mass density p, group velocity v group , and frequency intervals Ai/$. T(v) represent the 
transfer functions, i.e. the square of the transmission functions in Fig. [lb. Note that the 
transfer functions express the smoothing by radiative transfer for small wavelengths and we 
correct our measurements from this effect. From the GR and IGR models, we adopt the 
values pgr — 5.9 x 10~ 5 kgm~ 3 at 325 km and Pigr = 10 -4 kgm -3 at 175 km. We limit 
our study to the acoustic frequency range, i.e. from 5.2-15.1 mHz, the latter value being the 
Nyquist frequency of our observations. 

Wavelets are an appropriate tool to retrieve information on the distribution of power 
in space and time. With this technique we calculated 2D maps of power spectra in period 
bi ns of 20s from 70s up to 190s, i.e. the acoustic cutoff period. We employed the code 



by iTorrence fc Compel (119981 ) with Morlet wavelets and a level of significance of 95%. An 
analogous equation to Eq.[TJis then used in terms of periods U (with AUi =20 s) to determine 
the acoustic flux from the wavelet power averaged over time and FOV. 



Figure [2] shows averages of velocity flux spectra as functions of period (frequency) from 
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the wavelet and Fourier analyses. We differentiate between GR (thin) and IGR (thick) 
ffux before (dashed) and after (solid) correction of atmospheric transmission, in the case of 
exponentially increasing amplitudes. This correction is more conservative than that for waves 
with height-independent amplitude (see Fig.[Tb), therefore it will provide us with lower-limit 
corrected estimates. For the differentiation we obtain the power above GRs and IGRs from 
the wavelet analysis (cf. Fig.[3l). The resulting weights are transferred to the power in the 



Fourier analysis (see lBello Gonzalez et al.ll2010l ). 



Before correction, we find in both the Fourier and wavelet analyses, an IGR flux larger 
than the GR flux by a factor 1.7. The correction increases the flux in the IGR more strongly, 
as expected from the transmission functions (Fig. [lb). After the correction, the ratio between 
IGR and GR fluxes amounts to 1.8-2.0. 

TablaTJ 



The results on the total acoustic flux summed over all frequencies are collected in Ta- 
bled! The 'average' values refer to the flux over the FOV filled by GR over 46% of the area 
and by IGR over the remaining 54%. Results from the wavelet and Fourier analyses agree 
for the most part at low frequencies, i.e. for periods > 100 s. They amount to ~5500 Wm -2 
before correction and to ~ 6400-7700 Wm -2 after the conservative correction. Yet, a large 
difference appears at high frequencies, where the treatment of the Fourier and wavelet anal- 
ysis is different: (1) The wavelet approach is more restrictive in terms of noise when using 
a statistical significance level of 95%, therefore the retrieved values at these frequencies are 
low, (2) the wavelets were applied in period bins which reduces the spectral sampled points 
at high frequencies to just the bin at 80s and half of the bin at 100s, (3) the values at high 
frequencies are strongly corrected, i.e. the originally larger values in the Fourier analysis 
spectra are highly increased. 

We are confident in the correction of the flux spectra for frequencies < 10 mHz. The 
large correction of flux for frequencies > 10 mHz above IGRs from the Fourier analysis, see 



the in set in Fig. [21 turns the flux spectrum upward. This also was noted by iFleck et al. 



([20081 ) who point out the difficulty in determining the noise level. In addition, aliasing from 
power at frequencies above our Nyquist frequency may play a role. With these caveats, we 
consider the large fluxes at high frequencies as upper limits. The high frequency waves are 
most prone to radiative damping. Their amplitudes will thus be strongly reduced at the 
base of the chromosphere. 



In any case, even before correction, the directly measured energy fluxes are larger than 
those found in any previous observational study of acoustic waves. The nature of these 
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results and their consequences are discussed in Sect. HI 



3.2. Spatial distribution of the acoustic flux 



The wavelet analysis provides maps of wave power for each time step. With this tech- 
nique applied to the high-spatial resolution IMaX data, we can study in detail the distribu- 
tion of power (flux) over the observed FOV and its evolution with time. Figured shows an 
example of different levels of acoustic power over the period range [150-190] s overlaid on 
continuum images, and an animation on the time evolution of such power patches is available 
at URL. 

FigSJ 

The time lag corresponding approximately to the travel time from the bottom of the 
photosphere to 250 km height is about 42 s. The cadence of the time series is 33 s, therefore 
no delay between intensity and power maps has been taken into account. 

The wave power, far from being ho mogenously distr i buted , appears intermittently in 
spac e and time as already found by, e .g., iRimmele et al.l ( 119951 ). IWunnenberg et al.l (120021 ) 
and iBello Gonzalez et al.l (120091 120101 ) . Histograms of the acoustic power occurrence over 
intensity represented in Fig. [3] show that most of the power, at all frequencies, is associated 
with dark structures, i.e. structures with intensities I C ont/(^cont) < 1- We have referred to 
these structures as intergranules, but close inspection of the power maps overlaid on the 
continuum intensities (see animation available at URL, attached to this manuscript) reveal 
that significant power is found in forming granular (dark) dots and lanes. A second, less 
prominent component or shoulder is seen in the histograms at J C ont/ (leant) ~ 1.1. It has not 
been reported before, we see it thanks to the high contrast of the IMaX images. This power 
is found in small rapidly evolving granules, and in large granules (1) prior to their splitting 
and (2) at their bright borders. Yet not all splitting granules and bright borders show a 
wave power component. 



3.3. Energetic fast-fluctuating features 

Studying the evolution of acoustic power reveals patches of intense signal. They appear 
intermittently and randomly distributed, covering areas of typically 1 arcsec 2 and with power 
at all periods lasting for about 5 min. Such strong signals are found to be generated by rapid 
fluctuations in velocities up to ~4kms , that are of a diverse nature: (1) sudden (within 
2-3 min) up-down-up flows, (2) sudden down-up-down flows, (3) some associated to magnetic 
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field concentrations, (4) some with no significant magnetism, (5) some of the latter related 
to imploding granular areas, etc. 

In a first analysis, we have identified 32 of these events within the 23 min series. Masking 
these events and calculating the acoustic flux from their wave power, we find that they carry 
20-23 kWm -2 (estimated values before and after correction, respectively). Yet, their small 
area coverage, 0.1% of the FOV, implies that their global contribution to the total flux is 
not significant. 

Nevertheless, the identification of these energetic features at small scales from their 
strong power signature has opened a new field of study and we refer to a forthcoming 
contribution on the nature of these events. 



4. Discussion and conclusions 

In this study we have taken advantage of the high quality data of the IMaX/SuNRlSE 
instrument. The high spatial resolution provided new exciting findings, which (a) confirm 
that acoustic waves play an important role in the energy transport through the quiet solar 
atmosphere, (b) shed new light on the spatial and temporal distribution of propagating wave 
phenomena and (c) reveal the existence of very energetic processes occurring at small scales. 

Table[2] summarises the values of acoustic energy flux measured by different authors 
in the last years. Averaged values from the present study before and after correction of 
the transmission of the solar atmosphere to the velocity wave amplitudes are also listed for 
comparison. 

Tabled 

We note that the use of transmission functions for correcting velocity a mplitudes at 



short periods is controversial. The results from 3D hydrodynamic simulations by Fleck et al. 



(120101 ) suggest that the observed high frequency Doppler shifts are caused by rapid variations 
of the contribution functions in an atmosphere with high velocity gradients. These authors 
conclude that one should re-evaluate the cause of Doppler shifts with periods < 70 s, i.e. 
periods below the Nyquist period of our observations. 

We restrict the following discussion to the flux values from the wavelet analysis, which 
we consider to be lower limits. The average acoustic flux at 250 km is F aC)to t ~ 6 500 W m~ 2 , 
after a conservativ^l correction. It is mainly found at low frequencies (periods > 100 s). Yet, 



2 conservative stands for the less restrictive correction when calculating transfer functions for exponentially 
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the uncorrected flux already amounts to 5 500Wm -2 . T hese large values were never re- 
ported before, probably due to limited spatial resolution. IWedemeyer-Bohm et al.l (120071 ) 
pointed out from 3D simulations that a significant amount of acoustic power is lost when 
degrading the spatial resolution. This effect can be seen when appl ying to our data a 2x 2 
pixel binning. The signal is reduced by a factor of 1.4. We note that ICarlsson et al.l (120071 ). 
analysing HINODE data, found a lower p ower reduction upon binnin g. Also compare the 
resulting acoustic flux with that found by iBello Gonzalez et al.l (120091 ) at the same height. 
There the spatial resolution is at least a factor 2 less, although the spectral sampling is 
better. After transmission correction, the IMaX data provide fluxes larger by a factor of 2. 
We conclude that much of the power in the quiet Sun occurs at small scales and high spatial 
resolution is mandatory to detect it. 



Simulations by IStraus et al.l (120081 ) give an estimate of 30 % for energy losses between 
250 km up to 500 km due to radiative damping. Applying the same correction, the averaged 
acoustic flux we observe at 250 km would then be reduced to 4500 Wm -2 at the bottom of 
the chromosphere. On the one hand, the flux reduction between 250 km and 500 km may be 
larger because of the very small spatial scales in the present data set. On the other hand, 
the estimated flux at 250 km is only a lower-limit value since we have considered: (i) wavelet 
power, (ii) the conservative correction. This acoustic flux is comparable to the radiative 
losses of 4600 W m -2 qua ntified from the standa rd semi-empirical hydrostatic model of a quiet 
average chromosphere by lVernazza et al. (jl98lh . Chromospheric models calculated including 
radia tive losses from Fe II lines, arrive at energy needs of 14 000 W m -2 (lAnderson &: Athay 
19891 ). In our wavelet analysis, the contribution from high-frequency (>10mHz) waves is 
missing, but likely present in the solar atmosphere as suggested, e.g. the Fourier analysis. 

Let us recall the main findings of this first study on acoustic waves from IMaX/SuNRlSE 
observations at disc centre: 



1) The improvement in at least a factor of 2 in spatial resolution and the image quality 
of IMaX/SuNRlSE data have revealed a total acoustic power larger by at least a factor of 2 
than the largest power observed before. In view of the results obtained on the total acoustic 
power, we conclude that, at photospheric level, acoustic waves with periods >100s carry 
at least half of the flux nee ded to balance the o bserv ed radiative energy losses of the quiet 
chromosphere according to lAnderson fc Athayl (119891 ) . The possibility that the quiet solar 
chromosphere is supplied with acoustic wave energy to cover its radiative losses, as suggested 



increasing velocity amplitudes, Sect. 12.11 
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many years ago by L. Biermann and M. Schwarzschild and later on pursued intensively by 
P. Ulmschneider, needs to be revived and intensively studied, probing which fraction of the 
transported energy reaches the chromosphere. 



2) The main contribution to the acoustic energy flux is found in intergranules, as is 
extensively known. In addition, forming dark dots and lanes above splitting granules as well 
as small fast-evolving granules and the brighter borders of large granules are found to be 
sources of propagating waves. These are regions of str ong dynamics (turbu lence) as pointed 



out in the studies on the same IMaX data series by ISteiner et al.l (120101 ). and iRoth et al. 

d2oioh . 



3) Thanks to the high quality of the IMaX data, it has been possible to detect intermit- 
tent sources of strong energy fluxes at all frequencies in the acoustic domain. A preliminary 
analysis of their nature shows that they are generated by rapid (2-3 min) fluctuations in 
their LOS velocity occurring at small scales, i.e. typically of 0'.'3 x ff!3.. They show a diverse 
origin, in some cases they are related to small magnetic field concentrations, in others to 
imploding granules, etc. We refer to a forthcoming contribution on a detailed analysis of 
these events. 
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Fig. 1. — a) Response functions for velocity RF v (z) at line minimum of Fe I 525.02 nm: solid 
and dashed, for granule (GR) and intergranule (IGR) atmospheres, respectively; thick, after 
convolution of line profiles with spectrometer function and Gaussian fit, to mimic the IMaX 
observed profiles, b) Transmission of solar atmosphere to wave (velocity) amplitudes vs. 
wave frequency after convolution with the spectral transmission of the IMaX instrument; 
solid: GR, dashed IGR; thick: waves/ with amplitudes constant in height; thin: waves with 
exponentially increasing amplitude. Dotted lines indicate the frequency (acoustic) range un- 
der study, i.e. from the acoustic cut-off frequency, u ac =5.2 mHz (U=192 s) up to z/=14.3 mHz 
(U=70s). 
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Fig. 2. — Velocity flux spectra from wavelet (vs. period) and Fourier (vs. frequency) analysis. 
Thin and thick represent GR and IGR spectra, before (dashed) and after (solid) correction 
for atmospheric transmission, respectively. 



-14- 



Table 1: Total acoustic energy fluxes, F aCtto t [Wm 2 ] 



wavelet: 


uncorr 


corrected 


GR 


4430 


5000-4970 


IGR 


6400 


7725-8630 


average 


5495 


6470-6945 


5.2-10 mHz 


5455 


6360-6765 


10-14.3 mHz 


40 


110-180 


Fourier: 


uncorr 


corrected 


GR 


4930 


7000-6865 


IGR 


7540 


13075-20165 


average 


6340 


10280-14050 


5.2-10 mHz 


5630 


7720-8900 


10-14.3 mHz 


710 


2560-5150 



Note. - Left and right values in column corrected refer to the exponentially increasing and height- 
independent velocity amplitude approaches, respectively fSect. P?T2)) . 
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cont/ ^ cont^ 

Fig. 3. — a) Snapshot on acoustic power averaged over the period range [150-190] s, overlaid 
to a continuum image. The area corresponds to the upper-right quadrant of the full FOV; 
filled contours represent 12% (blue), 20% (green) and 32% (red) of the maximum power 
value over the sequence; tickmarks in arcsec. The underlying animation is available at URL. 
b) Occurrence of wavelet power of velocities above 1% of maximum power vs. normalised 
continuum intensities I CO nt/ {hont) for reconstructed data. Curves for [110-150] s and [70- 
110] s period bins were amplified by a factor 3 and 9, respectively, for better visualisation. 
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Table 2: Summary of observed total acoustic flux F aCitot [Wm 2 ] . 



method / data 


uncorr. 


corrected 


wavelet: binned (2x2) 


3975 


4600-4900 


wavelet: binned (lxl) 


5495 


6470-6945 


Fourier: binned (lxl) 


6340 


10280-14050 


other work 


Bello Gonzalez et al. (2010) 






h = 500 - 600 km, 0"4 




1700-2000 


Bello Gonzalez et al. (2009) 






h = 250 km, 0"4 


2500 


3000-3650 


Straus et al. (2008) 






h = 250 km, ~ 0"4 




1400 


h = 500 km, ~ C'5 




1000 


Carlsson et al (2007) 






h = 200 km, ~ 0'/22 




800 


Fossum & Carlsson (2006) 






h = 430 km, ~ 1" 




510 


Wunnenberg et al. (2002) 






h = 600 km, ~ 0'.'5 




900 



Note. - 



Notation in column corrected as in Table [T] 



